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Berry phases and the quantum-information theoretic notion of fidelity have been recently used to 
analyze quantum phase transitions from a geometrical perspective. In this paper we unify these two 
approaches showing that the underlying mechanism is the critical singular behavior of a complex 
tensor over the Hamiltonian parameter space. This is achieved by performing a scaling analysis of 
this quantum geometric tensor in the vicinity of the critical points. In this way most of the previous 
results are understood on general grounds and new ones are found. We show that criticality is 
not a suffcient condition to ensure superextensive divergence of the geometric tensor, and state the 
conditions under which this is possible. The validity of this analisys is further checked by exact 
diagonalisation of the spin-1/2 XXZ Heisenberg chain. 
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Introduction. - Phase transitions at zero-temperatures 
are dramatic changes in the ground-state (GS) properties 
of a system driven by quantum fluctuations. This phe- 
nomenon, known as quantum phase transition (QPT), is 
due to the interplay between different orderings associ- 
ated to competing terms in the system's Hamiltonian [lj]. 
Traditionally such a problem is addressed by resorting to 
notions like order-parameter and symmetry breaking i.e., 
the Landau-Ginzburg paradigm [2|] . In the last few years 
a big effort has been devoted to the analysis of QPTs from 
the perspective of Quantum Information Q the main tool 
being the study of different entanglement measures [J ■ 

More recently an approach to QPTs based on the 
quantum- information concept of fidelity has been put for- 
ward [f| , [|| . The strategy there is differential-geometric 
and information-theoretic in nature: GSs associated to 
infinitesimally close parameters are compared i.e., their 
overlap evaluated. The intuition behind is extremely sim- 
ple: at QPTs even the slightest move results in a major 
difference in some of the system's observables, in turn this 
latter has to show up in the degree of orthogonality i.e., fi- 
delity, between the corresponding GSs. Systems of quasi- 
free fermions have been analyzed 0,0 as well as QPTs 
in matrix-product states Finite-temperature exten- 
sions have been also considered showing the robustness of 
the approach against mixing with low excited states [Io| . 
Remarkably the fidelity analysis have been successfully 
carried over for the superfluid-insulator transition of the 
Hubbard model [HI]; this suggests that this framework, 
besides its conceptual appeal, can have some practical 
relevance even for fully-interacting systems where a sim- 
ple description is not possible. 

In Ref [13] it has been shown that the fidelity approach 
can be better understood in terms of a Riemannian met- 
ric tensor g defined over M. . Loosely speaking the singu- 
larities developed by g in the thermodynamic limit corre- 
spond to QPTs 12]. Even for finite-size systems, critical 
points have a markedly distinct (finite-size scaling) be- 



havior from the regular ones; in all the example studied 
so far this difference amounts to an enhanced orthogo- 
nalization ratio as a function of the system's size at the 
QPTs. Another intriguing relation between QPTs and 
geometrical objects i.e. Berry-phases, was suggested in 
Refs [HI and 14 1 . There it was argued that loops in the 
parameter space, encircling a critical line give rise to a 
non zero GS Berry phase even for an arbitrary small loop 
size. This fact indicates that at the critical points the 
curvature of the Berry connection should display some 
sort of singularity Q. 

In this paper we shall show that these two approaches 
share the same origin and can be therefore unified. We 
will perform a scaling analysis that allows one to un- 
derstand, from a single perspective, most of the results 
obtained so far in the fidelity approach and to investigate 
somewhat unexpected new ones. 

Geometric tensors. - We now lay down the formal set- 
ting. For each element A of the parameter manifold 
M. there is an associated quantum Hamiltonian H(X) = 
En= K ^«(A)|«'„(A))(* n (A)|, {E n+1 > E n ) acting over a 
finite-dimensional state-space TL; the mapping A — ► H (A) 
is assumed to be smooth. If | (A)) denotes the unique 
ground state (GS) of H(X) then one has the mapping 
^o'-M — > H: A — > |* (A)). More properly we will con- 
sider this map as a function valued in the projective space 
PTi. of rays. If the spectral gap of H (A) above the GS 
is bounded away from zero over M then Wo is smooth 
(l5| . The projective Hilbert space is the base manifold of 
a U(l) fiber bundle (l6| and it is equipped with a com- 
plex metric given by g(u,v) = (u, (1 — |^)(^ r |)u) (u and 
v denote tangent vectors in |^')(4 r |). Pulling this metric 
back to M. by \&o i.e., evaluating it on vectors of the form 
d/dt^>o(\(t)) one obtains the complex hermitean tensor 

Qvu := (d^old^o) - (0 M *o|*o><*o|$,*o) (1) 

Here the indices /i and v are labeling the coordinates of 
M. i.e., /i, v = 1, . . . , dim A4. This quantity is the quan- 
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turn geometric tensor (QGT) (l7j , both its real and imag- 
inary parts have a relevant physical meaning. 

The real part g^ v := ^RQ^ is a Riemannian (real) 
metric tensor over M. which defines the line element as 
ds 2 = ^^gfivdXfidXu- This metric has been shown to 
provide the leading term in the expansion of the fidelity 
between two GSs associated to slightly different Hamilto- 
nians [12]. More precisely if F(X, A') := |(#o(A), #o(A'))| 
is the fidelity then F(X, X+6X) w l-<5A 2 /2 g(d^ , d^ ) = 
1 — ds 2 /2 (i.e., g^v is the Hessian matrix of T[X 1 A') as 
a function of A' in A = A'). The meaning of this dis- 
tance function between parameters should be obvious: it 
is the Hilbert-space one between the corresponding GSs. 
This latter quantifies the operational distinguishability of 
two states [l8|; therefore even the induced metric g con- 



and the above inequality 



veys a definite information-theoretic meaning [12]. The 
fact that at the QPTs g exhibits singularities is consis- 
tent with the intuition that at critical points one has a 
major change in the GS structure i.e., it becomes "more 
different" , and makes it quantitative. Now we consider 
the imaginary part of {T]) F M „ := 3Q M „. Since the terms 
\P) are -from normalization- purely imaginary, one 
finds QQ^ = %{d^\d v V) = (0 M *|0„tf > - (0„* = 
dpAu-dvA^ where := {V\d v 9) is, for = |¥ (A)), 
the Berry adiabatic connection [191 ]. From this one sees 
that S(3 M „ is nothing but the curvature 2-form, respon- 
sible for the appearance of the Berry geometrical phase 
[3]. Of course for systems with real GS F is zero and 
the QGT coincides with its real part g. 

The QGT (pj can be cast in a way useful for later 
derivations as well as to decrypt its physical mean- 
ing. By inserting in Eq. ([T]) the spectral resolution 
1 = E„™o l*n(A))(*„(A)| and differentiating the eigen- 
value equation H(X)\^ Q (X)) = E (X)\^ (X)), one finds 
the identity 



y, (^o(A)|^g|^ n (A))(^n(A)|^g|W (A)) 

n=£0 



[E n (X) - So (A)] 2 

(2) 

This expression clearly suggests that at the critical 
points, where one of the e n (A c ) = E n (X c ) — Eo(X c ) > 
vanishes in the thermodynamic limit, the QGT might 
show a singular behavior. This heuristic argument is the 
same one proposed in ((l2j|) for the Riemannian tensor 



g^u and for the Berry Curvature F^ v in 14j . One of the 



aims of this paper is to establish this argument on more 
firm grounds. A similar scaling analysis in connection 
with local measures of entanglement at QPTs has been 
presented in (20j . 

We would like first to demonstrate an inequality useful 
to establish a connection between the tensor g and QPTs. 
We consider a system with size L d (d is the spatial dimen- 
sion). Since Q(X) is an hermitean non-negative matrix 



one has \Q 



< 



IQIloo = (4>\Q\4>) where \cj>) = {<t>X=x M 



denotes the eigenvector of Q corresponding to the largest 
eigenvalue. We set SH = ^^(d^H) then from Eq. (|2|) 



\Q,A < ^ £ - 2 |(* |^|*n)| 2 < £ i- 2 ^K*o|^|* n )| 2 
n>0 n>0 

= e ?((8H5rf)-\(8H)\*), (3) 



where the angular brackets denote the average over 
|\&o(A)). Now we assume that the operator 5H is a 
local one i.e., SH = £^,-<5V^; then the last term in 

Eq. © reads E^W/) - (6VJ(6V?)) If the GS 
is translationally invariant this last quantity can be writ- 
ten as L d J2 r K ( r ) ■= LdR where K ( r ) ( 6V * SV i+r) ~ 
(5Vi)(SV^ +r ) independent on i. For gapped systems i.e., 
£i(oo) := limi^oo £i(L) > the correlation function 
G(r) is rapidly decaying [2lj and therefore K is finite 
and independent on the system size. Using ([3]) it follows 
that for these non-critical systems \Qn V \ cannot grow, 
as a function of L, more than extensively. Indeed one 
one has that lim^^oo \Q^\/L d < Ke^ 2 ((X)) < oo. Con- 
versely if liniL^oo \Q^\/L d = oo i.e., \Q^ V \ grows super- 
extensively, then either £\{L) — > or K cannot be finite. 
In both cases the system has to be gapless [22j. Sum- 
marizing: a super-extensive behavior of any of the com- 
ponents of Q for systems with local interaction implies a 
vanishing gap in the thermodynamic limit. 

This sort of behavior has been observed in all the sys- 
tems analyzed in Refs and does amount to the 
critical fidelity drop at the QPTs. As we will show 
in the next section, in general the converse result i.e., 
QPT^ super-extensive growth of ds 2 (L) does not hold 
true: Q^/L d can be finite in the thermodynamic limit 
even for gapless systems. In order to demonstrate this 
fact we now move to a representation of in terms of 
suitable correlation functions. This key move is an ex- 
tension of the results of You et al. [23| , for the so-called 
fidelity susceptibility, to the whole QGT. 

Correlation Function representation. - Let us consider 
the following imaginary time (connected) correlation 
functions 

GV(r) = 8(r) ({dpH(T)duH(0)) - (d„H (0))(d„H(0))) , 

(4) 

where X(t) := e T Xe rH . Using again the spectral 
resolution of the identity associated to -ff(A) one finds 
G^ir) = e(T)J2 n>0 e- e ^X nfl X* v , where X nfl := 
(* (A)|9 M -ff|*„(A)). Notice that if H(X) = H + XV 
then G(t) is nothing but the dynamic response func- 
tion associated to the "perturbation" V. We now move to 
the frequency domain G M „(o;) = /_ °° dre~ lU17 G m „(t) = 
J2 n >o X rifiX*^(iuj + £„) _1 By comparing this equation 
with © it is immediate to see that 



d ~ f +oc 
-i— Guv(u)\u=o = / drrGuu(r) 



(5) 



This equation is an integral representation of the QGT in 
terms of the (imaginary time) correlation functions of the 
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operators d^H. Eq ((4]) is remarkable in that it connect the 
tensors g M „ and F^ v directly (and non perturbatively) to 
the dynamical response of the system to the interactions 
<9 M -ff's. In this way geometrical and information-theoretic 
objects F and g are expressed in terms of standard quan- 
tities in response theory and their physics content is so 
further clarified. Eqs. (|4]) and (0 provide the starting 
point for our scaling analysis 

Scaling behavior.- First we assume that the operators 
d^H are local ones i.e., d^H = J2 X V^(x).We also rescale 
the QGT (TJ) by the system size Q^ v — > q^ = L~ d Q llv 
in order to obtain well-defined quantities in the ther- 
modynamic limit. Now we consider the scaling trans- 
formations x — > ax, t — > ofr, (a E R + ). Assuming 
that, in the vicinity of the critical point A c , the oper- 
ators Vfi have well-defined scaling dimensions [24| one 
has V M — > a "Vp, these relations along with Eqs. (|4]) 
and ||5j) imply 



A Q — A 



2(-d. (6) 



For simplicity we assume now that there is only one driv- 
ing parameter A and drop the indices ijl and v in A^„. 

If £ is the correlation length one has £ = | A — A c | ^ 
(here v is the correlation length critical exponent) and, 
if A>, is the scaling dimension of the driving parameter 
A, v = A^ 1 . Putting all this together and following stan- 
dard arguments in scaling analysis one obtains that (in 
the off-critical region, L ^> £) the singular part of the 
-intensive- QGT behaves as 



<7^(A « A c ) ~ |A - A c | 



a q /a> 



(7) 



Instead at the critical point i.e. £ = oo, where the only 
length scale is provided by the system size itself one gets 



?/«/(A = A c ) 



(8) 



Equations ([7|) and ((HJ represent the main result of this pa- 
per. From Eq. (|7|) one sees that close to the critical point 
the QGT g M y is diverging for Aq/A\ < 0; on the other 
hand when one is sitting exactly at the critical point one 
finds that, besides an extensive contribution coming from 
the regular part, the singular part contributes to Q^ v in 
a manner which is: i) super-extensive if Aq < ii) exten- 
sive if Aq = and iii) sub-extensive for Aq > 0. Hence 
we observe that q^ can be finite at the critical point, 
even in a gapless system, provided Aq > 0. An explicit 
example of this phenomenon will be discussed in the se- 
quel; before doing that we show that this analysis allows 
one to understand in a unified manner the results found 
for quasi- free fermionic models in ln quasi- free 

fcrmionic models the most relevant operator admissible 
has scaling dimension equal to one, therefore from (|7|) 
and ((8]) one finds, close to A c , Q M „ = g^ v ~ L\ \ — A c | _1 
and Q M „ ~ L 2 at the QPT. Notice that if AQ < -1 one 
expects super-quadratic behavior. 



XXZ chain - We provide now a further test for the 
ideas presented in this paper: the S = 1/2 XXZ 
Heisenberg chain. The model is defined by H = 
JT,i [SfjSLi + SfSf+i + Afif 5f +1 ] . It is well known 
(see e.g. [25j) that in the regime A 6 (—1,1) the model is 
in the universality class of a c = 1 conformal field theory 
(rf = C = 1) displaying gapless excitations and power low 
correlations. For A > 1 the model enters a phase with 
Ising-like antiferromagnetic order and a non-zero gap. Fi- 
nally, the isotropic point A = 1 is a Berezinskii-Kosterlitz- 
Thouless transition point. The low energy effective con- 
tinuum theory is given by the sine-Gordon model: H 



(yi6^K^ here 



is the bare Fermi velocity, u the renormalized one, K is 
related to the compactification radius of the field (f> and 
a is the lattice spacing (we use the notations of [HI)- 

We now analyze the behavior of the fidelity when the 
anisotropy parameter A is varied. Correspondingly we 
are interested in the operator V (x) — S^S^. +1 and its 
scaling exponents. In the continuum limit the operator 
V (x) contributes with a marginal operator - (d x $>) 2 - 
with scaling exponent Ay = 2, plus a term which is 
precisely the cosine in the sine-Gordon Hamiltonian. The 
cosine term is irrelevant for |A| < 1, marginal at A = 1 
and relevant for A > 1 where it is responsible for the 
opening of a mass gap. Its scaling dimension is precisely 
equal to AK. The parameter K can be fixed by the long 
distance asymptotic of the correlation functions obtained 
by Bethe Ansatz, and one gets, for A 6 (—1, 1] K = 
ir/2(ir — arccos (A)) -1 . 

When one considers the finite size scaling of q in the 
gapless region |A| < 1, the leading contribution is dic- 
tated by the most relevant component of V (x) which in 
this case is the marginal one i.e. Ay — 2. Correspond- 
ingly, using equation 1(5)) we obtain Aq = 1, so that, in 
the whole region |A| < 1, the finite size dependence of 
the (rescaled) QGT tensor is q = Ai + A^L -1 , where 
A\p are constants which depend only on A. Note that 
the term A\ is the contribution coming from the regular 
part. This kind of scaling has to be contrasted with the 
one observed in quasi-free fermionic systems, where one 
has q = A' i + A'%L. We would like to stress again that 
in general one expects a super- extensive behavior of Q, 
and a corresponding fidelity drop when d + 2£— 2 Ay > 
i.e., when the operator associated to the varying param- 
eter is sufficiently relevant. When this condition is not 
fulfilled the rescaled QGT tensor does not diverge in the 
thermodynamic limit at critical points. Nevertheless a 
proper finite size scaling analysis allows one to identify 
the critical region. To check this latter feature as well as 
the predicted scaling behavior we have performed exact 
Lanczos diagonalizations. The agreement between nu- 
merical data and the theoretical prediction is shown in 
Fig ©. 
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g/L o.oi; 




its quantum-critical behavior. 
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Figure 1: Scaling behavior of metric g. The points are the 
data obtained via Lanczos diagonalization of periodic chains 
of length L up to 26, using equation ^(A, A + <5A) « 1 — gS\ 2 /2 
with 5\ = 1.0 x 10" 3 . The points at |A| < 1 are well fitted 

(solid lines) with g/L = Ai + A 2 L~ Y + A 3 L 3 ~ 2A v . The con- 
tribution A% comes from the irrelevant operator with scaling 
dimension Ay' = 4K. As this operator becomes rapidly irrel- 
evant for |A| < 1 its contribution can be hardly observed. At 
A = 1 a better a fit is obtained with logarithmic corrections 
as expected at the isotropic point. In the massive regime we 
expect that the thermodynamic limit is approached exponen- 
tially fast. We obtain a good agreement by fitting our data 
with the phenomenological formula g/L = A\+A2e~ L ^L~ 1 ^ 2 
where £ is the correlation length as given by Bethe Ansatz 
As £ (A = 2) = 8.35 ... we used only points with L > 14. 



Conclusions.- The mapping between a quantum 
Hamiltonian and the corresponding ground state endows 
the parameter manifold with a complex tensor Q. The 
real part of Q is a Riemannian metric g while the imag- 
inary part is the curvature form giving rise to a Berry 
phase. The metric tensor g is closely related to the quan- 
tum fidelity between different ground states; in the ther- 
modynamical limit it has been shown to be singular at 
the critical point for several models featuring quantum 
phase transitions e.g., quasi-free fermionic models. The 
same kind of singularity have been argued to exist for the 
form F and the associated Berry phases as well. In this 
paper we demonstrated that i) the components of Q 
have an integral representation in terms of response func- 
tions ii) a super-extensive behavior of any of the Q^'s 
implies, for local models, gaplessness iii) the singular part 
of the Q^y fulfills scaling relations which explicitly con- 
nect their singular behavior with the universality class 
of the transition i.e., critical exponents. In particular 
these relations show that gaplessness is a necessary but 
not sufficient condition for a super-extensive scaling of 
the metric tensor, i.e. enhanced orthogonalization rate. 
The theoretical analysis has been further supported by 
a numerical study of the finite-size scaling of the fidelity 
for the XX Z spin 1/2 chain. The main message of this 
paper is that apparently unrelated results can be un- 
derstood in unified fashion by unveiling the underlying 
common differential-geometric structure and analyzing 
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